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As is known, the free heat-kernel on the integers (a modified Bessel 
function) is turned into the periodic free heat-kernel on the discrete 
circle by factoring, giving a pre-image sum. I generalise existing treat- 
ments by making the functions periodic up to a phase, thus introducing 
an extra parameter into the analysis. Identifying the classical paths 
form with the conventional eigenfunction expression, I find a com- 
binatorial trace identity which allows various Bessel identities to be 
extracted, such as a generalisation of the Jacobi-Anger expansion. 
The free Dirichlet, Neumann and hybrid Dirichlet-Neumann heat- 
kernels on a discrete interval are constructed using both modes and 
images. The Neumann imaging mirror has to be placed at a half- 
integer. 

The corresponding lattice Green functions are expressed in terms of 
Chebyshev polynomials and the Laplacian matrices extracted. The 
generating functions for circuits with bumps are evaluated. 
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1. Introduction 



In the next two sections I summarise and rework some known material, [1] , [2] , 
on the free heat-kernels on the integers and on the discrete circle with p vertices 
(or p-cycle). In section 3, I extend the discussion to complex functions on the 
circle periodic up to a phase. This is mathematically possible and gives extra 
room for manoeuvre in the analysis and allows for a more general trace identity. 
Other sections contain some algebraic consequences, in particular Bessel function 
identities, which might be novel. I also construct the Dirichlet, Neumann and 
Dirichlet-Neumann kernels on a discrete interval using both modes and images and 
then, by Laplace transformation, Chebyshev polynomial forms for the lattice Green 
functions. 

I am only concerned here with one dimension but higher dimensions can be 
reached by cartesian products and, as physical motivation for revisiting this topic, I 
remark that the spectral problem, e.g. determinants, on higher dimensional lattices, 
particularly two dimensional ones, has relevance in statistical mechanics, [3], [4], 
string theory, [5], [6], network theory and many other areas. 

The two-dimensional discrete (periodic) lattice also appears in the description 
of the generators of SU(N), [7], [8], and can have the interpretation of a phase 
space. It occurs thus in a finite model of two-dimensional ideal hydrodynamics on 
the torus, [9], [10]. In [9], the discrete (^-function appears when constructing the 
stream function of a lattice vortex by analogy with the continuum version. 

2. Heat— Kernel on the integers 

The free heat-kernel, K^(j,j';t), on the integers, j,f G Z, has been given by 
Chung and Yau, [1], for example. The defining equation is the discrete heat, or 
diffusion, equation 2 




(1) 



with initial condition, 




for t < 0. 



(2) 



VAy(j) = y(j + 1) - 2y(j) + y(j - 1) . 



(3) 



Some authors have a factor of 1/2 in the first term. 
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The explicit expression for K is derived in [1] and [2]. I give an algebraically 
slightly different, and more rapid, development. 

By translation invariance, it is sufficient to set j' = and consider K z (j;t) = 
K z (j, 0;t). It is also expedient to remove the middle term in the Laplacian, (3), at 
the outset by setting, 

K z (j;t) = e- 2t K z V;t), 
so that the operator solution of (I) is effectively expressed in, 

K z (j; t) = e*^" 1 ) 8(j) = e tE 5{j) , (4) 

where E is the usual stepping operator, Ey(j) = y(j + I), and S(j) = <L , in terms 
of the Kronecker delta. 

Expanding the exponentials in (4) it is immediately seen that, 

« ur '*(>) = j, 

K Z (j;t) = e tEt - = J2 , T -„ J>0, 

which is the standard series expression for the Bessel function, 7 (2t), and I obtain 
the published formula for the free heat-kernel on (all) the integers, 3 

K z (j;t) = e- 2t I j (2t), t>0. 

This result is not surprising, cf [2], in view of the fact that Basset's modified Bessel 
function, I u , satisfies the recursion, 

(E + E- 1 )I u (z) = -^-I iy (z). 

Further, the initial condition, (2), is validated using i»-(0) = 8(j). 
By discrete translation invariance, 

K z VJ , ;t) = e- 2t I j _ j ,(2t). (5) 

Karlsson and Neuhauser [2], note that a probabilistic derivation of this result is 
given by Feller, [11], and also remark on the relevance of the composition rule, 4 

oo 

Ij-j"( Z l + Z 2) = Ij-j'( Z l) Ij'-j"( Z 2) ■ 

j'=-oo 

"^Equivalently, one requires the term E~ J in (E + E^ 1 ) n in order to pull the index j in Sj back 
to zero, cf [1]. 

4 By an infinite folding, one could obtain an analogue of a path integral. 
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As is very well known in the continuum case, the heat-kernel on the real line, 
K-^, can be found in many ways, the standard one being eigenfunction expansion. 
Equating these various expressions gives identities. I will not do this for the integers 
because I wish to pass immediately to the factor, Z/pZ, which is the discrete circle 
with p vertices (a p-cycle, C). obtained from the covering, Z, by identification, 
exactly as for the continuum. 



3. Heat— Kernel on the discrete circle 



A standard, geometric, way of obtaining the heat-kernel on the covered space 
is by a pre-image sum of that on the covering space to give a periodic function. 5 
Hence, 

oo 

K Z/ P zU> 3'\ t) = E f + mp;t) 

(6) 

Ij-j'+mp(2t) , 



m= — oo 

oo 

-It 



= e 

m= — oo 



which can be called the sum over classical paths form, by analogy with the more 
common continuum situation. I do not distinguish, notationally, between points on 
Z and those on Z/pZ. Quantities on the covering space are the thermodynamic 
limits of those on the covered space and correspond to the direct term, m = 0, in 
the sum (6). 

A dual representation arises from the eigenfunction expression. The Laplacian 
eigenproblem on the discrete circle is old, e.g. [12]. I therefore just write down the 
eigenfunctions, in complex form, and the eigenvalues. 

The normalised modes on the discrete circle are, 

Vn0') = ^e 2 ^>, n = 0,...,p-l, (7) 
VP 

with corresponding eigenvalues, 

A n =4sin 2 — , (8) 

so tha 

1 p_1 

^z/pzO';*) = -^e- 4sin2 (™/^V^'/P. (9) 

P n=0 



5 This procedure has occurred in numerous places at various times. 
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Equating the two representations of the heat-kernel yields the identity, [2] , 

oo 1 p— 1 

e~ 2t I j+m P ( 2t ) = ~ J2e~* sin2{7rn/P)t e 27vinj/P - (10) 

m= — oo ^ n=0 

Cosmetically re-expressed, this simplifies to, 

oo p-1 

E I j+mp (z) = ~Y1 e cos(2 ™ /p) z e 27rtnj/p . (11) 

m= — oo ^ n=0 

For j = 0, this result is to be found in Al-Jarrah, Dempsey and Glasser, [13], 
together with numerous Bessel series. Their method involves a generalised Poisson 
formula due to Titchmarsh. 

Since the left-hand side of (10) is real (if 2t is) one can conclude, 

oo 1 p-1 

e " 2t E W( 2 *) = -E e ~ 4Sin2(WP)t ^s2nnj/p 

m= — oo ^ n=0 fl2l 

= ^ e -4sin 2 (Wp)* S in27rnj/p. 

n=0 



4. Heat— kernel on the circle for twisted fields. 

So far, I have reproduced published results. I wish now to generalise the analysis 
a little. In the continuum case, according to the results of Schulman, [14], Laidlaw 
and DeWitt, [15], and Dowker, [16], in the quantum mechanical context the wave 
functions on the covered space can be be twisted by a unitary representation of 
the character group of the fundamental group, 7r l5 of the covered manifold. 6 For 
the circle, ir 1 = Z, and the unitary representations are generated by the phase, 
exp27ria, where a is a real parameter. The wavefunction picks up this phase as 
its angle argument increases by 27r, i.e. as the circle is completed. A possible 
interpretation of this involves an Aharonov-Bohm flux. 

In terms of the propagator, this twisting translates into the statement that 
each term of the pre-image sum enters with a multiplying factor of a power of this 
phase. 

6 For line bundles, tt 1 can be replaced by the more tractable first homology group, H 1 as discussed 
in [17], [18]. 



4 



Here I want to use this flexibility in the discrete case, and for diffusion. That is, 
I consider a (now complex) function, defined on the points, j, and satisfying 

the twisted periodicity condition, 

iP(p) = e 2 ™V(0) , ^(P - 1) = e 2 ™V(-l) , (13) 

which generalises the usual one. 7 At the moment, I present this simply as a math- 
ematical generalisation. 8 

Firstly, the classical paths expression is modified to the twisted form, 



oo 



m= — oo 

oo 



m= — oo 



Then, instead of (7) and (8) the eigenfunctions are, 

C0') = ^e 2 ^+^>, n = 0,...,p-l, 0<a<l, 



with corresponding eigenvalues, 



TV 



(n + a) 



X n (a) = 4 sin 2 — — - — '- . (14) 

As before, identifying the paths and eigenfunction expressions yields, instead 
of (10), the twisted identity, (z = 2t), 

oo 1 p-l 

e -z e" 2 ™ I j+mp (z) = ~J2 e" 22Sin2 l^ n+ ^/P) e 2*i(.n+a)j/p _ ( 15 ) 

m=— oo P n=0 

I record the special case a = 1/2, (anti-periodic functions), 

oo 1 p— 1 



]T (-1)™ ^+ mp (*) = - E e" 22Sin2 fa 2 ^ 1 )/ 2 *) e -(2n+lb/p 



m= — oo n=0 

p-l 

P 



1 £ e -2 2 sin2 ( 7r (2n+l)/2p) ^ ^ + 1)j/p) ^ 



n=0 

(16) 



used later. 

When j = see also the computations of Cojocaru, [20]. 



7 Fort, [12]. discusses the case a = 1/2 which corresponds to (real) anti-periodic functions. 

8 The use of such twisted boundary conditions is commonplace in statistical physics and string 
theory, both continuous and discrete. As a typical example I cite [19]. 
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5. Further identities 



Before discussing the above results, I outline some known things. The textbook 
identity, 

| | I 4 sin 2 (n(m + a) /p) + 4 sinh 2 7 J = 2 ( cosh 2p"f — cos 27ra) , (17) 
m=o ^ ' 

can be proved by elementary trigonometry, e.g. Hobson, [21], Loney, [22], §368, 

Bromwich, [23], p. 211, and so can reasonably be considered as given. 

It can be understood as equating two expressions for the determinant of the 
operator — VA + 4 sinh 2 7. The left-hand side is the product of the eigenvalues 
(the definition) and the right-hand side results from a non-eigenvalue evaluation 
using, say, the Gel'fand-Yaglom technique. This could be considered an alternative 
derivation to the trigonometric one. 

Taking logs and differentiating with respect to 7 produces the known, 

' ' 1 1 sinh 2^7 1 



_ y t . 

4f» ^ Q sin 2 (yr(m + a) /p) + sinh 2 7 2 sinh 27 cosh 2^7 - cos 2na 

t/ p _ 1 (cosh2 7 ) 



(18) 



2 (T p (cosh 27) - cos27ra) ' 

where T and U are Chebyshev polynomials of the first and second kind. 

In [24] I gave a fragmentary history and some discussion of these expressions in 
connection with Verlinde's formula for the dimensions of vector bundles on moduli 
spaces. In particular it was noted that the values of the twisted (^-function, 

P-i 1 

Cz/ P z( s ; «) = Yl (a ■ 2 ; 7 w ^ ' ( 19 ) 

11 (4 sin TT{m + a)/p) 

on the discrete circle, at positive integers, i.e. , 

p-i 

2 ~ 2n cosec 2n (n(m + a)/p) , a>0, (20) 

m=0 

can be obtained by expanding (18) in powers of 4 sinh 2 7. This is just an example of 
the Euler-Rayleigh technique. See also Dikii, [25]. 

The quantity 4 sinh 7 can be thought of as a mass squared, m , a Laplace 
transform parameter, s, a fugacity, z, or as a spectral parameter, A. One interpre- 
tation of the left-hand side of (18) is a resolvent which is the derivative of logdet 
with respect to m 2 . 
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Alternatively, setting 7 = in (17) and taking logs reduces to the ancient, and 
directly derivable, Kubert identity, 

p-i 

log 2 sin Tv(m + a)/p = log 2 sin ira , 

m=0 

which can be differentiated, this time with respect to a, also to give the cosec sums 
(20). For general interest, I give two examples, e.g. Fisher, [26], 

Cz/pZ(l' a ) = 4 ^cosec 2 7ra/p + p 2 cosec 2 a — cosec 2 a/p*j 

Cz/ P z( 2 ' a ) = Yg ^ cosec 4/Tra /P + f 4 cosec 4 a — -p 2 (p 2 — 1) cosec 2 a — cosec 4 a/p^ . 

The first term is to be dropped if a = as it corresponds to a zero mode and is not 
included in the ^-function. 

Other relevant formulae are the Fourier series, 

sin27ra 9 _ 

2 2_,e~ 27n sin27rna 



cosh 27 — cos 2-rra 

n=1 00 ( 21 ) 

- Slnh27 = l + 2y e - 2 ^cos27rna, 7 >0, 

cosh 27 — cos 2na ^-^ 

n=l 

again derivable by school trigonometry, e.g. Loney, [22] §§293-294. 

After this recapitulation, I can now return to the results of the earlier sections. 
Karlsson and Neuhauser, [2], derive the identity, (18), for a = 0, by taking the 
Laplace transform of the heat-kernel equation, (10), evaluated at j = (which 
amounts to taking the trace 9 ). The right-hand side gives the resolvent, correspond- 
ing to the left-hand side of (18). The transform of the summand in the left-hand 
side of (10) can be given as a standard closed form and, on summation, effectively 
gives the right-hand side of (18). 10 

It is no harder to show that the more general identity, (18), follows on Laplace 
transforming the heat-kernel relation, (15), at j = 0. As used in [2], in a different 
notation, the tabulated Laplace transform, 

POO -I 

J dze- z I n (z)e- 2zsinh2 "t = sinh27 e ~ 27n > n >°> ( 22 ) 



9 The factor of p is the volume 
10 The relation between variables is that the s of [2] equals i/2sinh7 here. 
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allows the summation over m to be performed immediately, after this has been 
written over non-negative integers using I_ n (z) = I n (z). Simple algebra then 
quickly delivers (18). Depending on one's starting point, this can viewed either as 
a derivation or as a confirmation of (18). 

Continuing with this theme, starting with the basic integral form of the modi- 
fied Bessel, 

I ( z ) = - I #e 2COS< ^cosn</>, 
n Jo 

the Laplace transform, (22), can be derived purely trigonometrically on applying 
the elementary Fourier series (21). Alternative proofs can be more complicated 
and may involve contour manipulation, or series expansion. (Consult Gray and 
Mathews, [27] Chap.VI and Ex.l4,p.76.) 

I finally note that the heat-kernel relation (15) can be more neatly written as 
in (11), 

00 p-l 

e -2irimaj_^ ^ _ ±_ e zcos (2Tr(n+a)/p) e 2-ni{n+a)j /p ^3) 

m= — oo ^ n=0 

Some consequences of this are derived in the next section. 

6. Bessel relations 

Taking the trace of (23) (i.e. setting j = 0) gives the Fourier cosine series, 

1 oo p-l 

-I (z) + ™s2vmal mp {z) = - £ e 2COS (M»+«)/p) . (24 ) 

m=l n=0 

(See also the analysis of Cojocaru [20].) 

For example, for p = 1 and p = 2, I find, 



1 1 

-I (z) + cos2nmaI m (z) = - e z cos2 ™ 

m=l 

1 oo 1 

-I (z) + cos 2imia I 2rn (z) = — cosh(z cos ira) . 

m=l 

Next, setting a = 1/2 in these produces the anti-periodic expressions, 



(25) 



-z 



^0W+E(- 1 ) m/ mW = ^- 

(26) 



1 00 1 

«W+E(- 1 ) mj ^)=9- 



m=l 



8 



Adding the first of this last set to the p = 1, a = equation in (25), gives the 
p = 2, a = relation in (25). Subtracting, on the other hand, yields, 

oo 

2 5>2 M+ i(*) = sinh*. (27) 

These relations can be obtained from those in [13], which are actually expressed 
in terms of the J n Bessel function. They are listed results, [13] providing the 
references. 

Furthermore, the first equation in (25) is recognised as a real form of the 
Jacobi-Anger expansion which, as is well known, yields a generating function for 
I m if one puts t = e 2lTlOL when it turns into the basic, 

oo 
m= — oo 

The other interesting value is a = 1/4 for which (25) becomes, 11 

oo 

i = / ( Z )+25](-iru) 

m=l 

oo 

cosh(^/v / 2) = I (z) + 2 £ (-l) m • 



(28) 



m=l 



The first of these agrees with the second of (26) . 

For comparison, I present an elementary derivation of the first equation in (28). 
First I note by simple cancellation and recursion that, (all arguments are z), 



h= E(- 1 ) m+1 Wl+E(- 1 ) m/ : 

m=0 m=l 

oo 

= E(-i) m (/ 2ra -i+w) 

m=l 

oo 

= E(- 1 ) m/ 2 



L 2m 

m=l 



Then, using I 1 = Iq/2, integration yields the required relation after setting z = 
to fix the unknown constant. 



11 The value a = 1/4 is a turning point between the 'boson' value, a = 0, and the 'fermion', 
a = 1/2. 
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Finally in this section, I give two further examples of (24), this time for p = 3, 
1 00 1 

2^0 (z) + E(- 1 ) m/ 3-^ = g(e 2/2 + 2e-*/ 4 cosh 3z/4) 

m=l 

1 °° 1 

2 + E W*) = 6 ( X + 2 C ° sh ^ /2 ) • 

m=l 

For other uses of similar Bessel relations see [28] . 



7. Heat— Kernel on the discrete Dirichlet interval 

It is straightforward to find the Dirichlet (D) heat-kernel on the discrete in- 
terval, or path, V, defined to be a set of p + 1 vertices, labelled by j G Z. The end 
points, j = and j = p, are boundary vertices at which the function, and hence 
the heat-kernel, vanishes. The remaining, internal vertices, j = 1, . . . ,p— 1, are the 
dynamical ones. Then V is a (p — l)-path. 

Translation invariance is lost and one has to consider K®(j,j';t) which, just 
as in the continuum case, can be constructed from K^(j, j';t) by images. The 
positioning of the continuum images is classic and given e.g. in Thomson, [29], and 
has been used frequently in heat conduction and hydro dynamical problems. (See 
Hicks, [30], Carslaw and Jaeger, [31], Basset, [32] §57.) In the discrete situation, 
the mirrors are also placed at j = and j = p, enclosing the p — 1 dynamical, free 
vertices and it is easy to see that the discrete Dirichlet conditions hold at the end 
points. 

This construction results in the image summation, 

oo 

K°(j,j';t)= {K Z (j,j'-2mp;t)-K z (j,-j'-2mp;t)) 

(29) 

= e E (lj-j'+2m P C2t) — Ij + j'+2mp(2t)) • 



m= — oo 



Dual to this image form is the eigenfunction expression. The eigenvalues and 
modes here are genuinely historic, being, 



2 ym 



A n = 4sin 2 — , n = l,...,p-l, (30) 
2p 



and 

/.x /2 . jirn 
VnW = \ Z sm 



p p 
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respectively, e.g. Fort, [12], so that the heat-kernel is, 



K°(j,j';t) = - Y e - 4t - 2 ("/2 P ) sin ^ sin !^ 

( 3i ) 

= i y e -4tsin"(W2p) ( cos ^(J-/) _ cos ™V + j') 

which can be seen to be equivalent to the image form, (29), by virtue of the identity, 
(12), given earlier, after setting p — > 2p. I give a few details, as they are not quite 
obvious. I first note that the n = term in (12) cancels by the subtraction in (31). 
Then I remark on the identity, 

p-i 

^cos(7mj/p)(e- 4tsin2 ^ n / 2p) - (_i)^-4tco S 2 (W2 P )) 

n=l 

= e- 2tS ^cos{nnj/p)(e 2tcos ™ /p - (_i)i e -a««Wpj =0 . 

n=l ^ ' 

which comes into play when the sum Yl^!Li m (12) is split into one from n = 1 to 
p — 1 plus one from p to 2p — 1 and shows that the two sums are equal so that (12) 
becomes, 

e~ 2t E ^ + 2 mp (2t) = - + -E e " 4Sm2(Wp)tcos2 ^'/P- (32) 

m= — oo n=l 

The required equivalence of (31) and (29) then follows easily. 

In the continuum case, this equivalence is effectively Poisson summation, as 
has been known for about 150 years, in various guises, and the identity (10) is a 
discrete analogue, [1], [2]. 

Taking the trace, i.e. setting j = j' and summing over j from 1 to p — 1, yields 
the image form, 

p — 1 oo oo p — 1 

j2 K vUJ-,t) = e- 2t (p-i) w 2 *)- e ~ 2t E E^W 2 *) 

j—1 m= — oo m= — oo j=l 

oo oo p— 1 

= e- 2t P E h mP m-e- 2t E J 2,W 2t ) 

m= — oo m= — oo j=0 (33) 

oo oo 

= e- 2t P w 2 *)- e ~ 2t 



m= — oo 
oo 



= e~ 2t p J 2m P (2t) - e~ 2t cosh2t 



m= — oo 
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where I have used the listed (25), for a = and have split all the integers into 
residue classes mod p, N = j + rap, m G Z and < j < p. 

Equating this to the trace obtained from the eigenfunction form yields an iden- 
tity that contains nothing new over the periodic the equivalence, proved 
just above, shows. 



8. Heat— Kernel on the discrete Neumann interval 

Again I approach the construction from the image side and now offset the 
mirrors by 1/2 to the half-integer points 1/2 and p + 1/2 so as, this time, to enclose 
a p-path, V (j = 1 to j = p) of p free vertices. The end points, j = and 
j = p + 1, which are sometimes referred to as fictional, or ghost, points, are now to 
be regarded as images and tied to the free vertices at j = 1 and j = p, respectively, 
thus reproducing the standard discrete Neumann conditions. In the continuum limit 
the 1/2 offset becomes irrelevant. 

It is easily concluded that the even images of the point j are at j — 2mp and 
the odd ones at —j — 2mp + 1 (m = — oo — > oo), so as to yield the image form, 

oo 

K»(j,j';t)= {K Z (j,f-2mp;t)+K z (j,-f + l-2mp ] t)) 

(34) 

= (K Z (j-j' + 2mp;t) + K z (j+j'-l + 2mp;t)), 

m= — oo 

both sets of images coming in with the same sign. 

I now look at the eigenfunction expression. The eigenvalues are as in (30) 
except that the mode label extends to n = 0, corresponding to a zero mode. 

The iV-eigenfunctions on the p-path are again standard, 



N ... 2 mr(2j-l) 

= Vp COS 2p ' "= 1 »-"P- 1 
1 

= — , n = 0, 
VP 

and the free, off-diagonal heat-kernel can be written, 

K?(j, j>- t ) = - + -Y e"« ^ (W2P) cos ™W ~ cos ™W ~ !) 

p p ^-^ 2p 2p 

n=l 
p-1 



(35) 



= - + -y e-^^(cos + cos . 

p p^ x v p p J 

(36) 
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To transform this into an image form one proceeds as in the last section by 
using (32) to obtain, 

oo 

K^(j,j';t) = e 21 ^2 {lj-j'+2m P (2t) + ij + j/_i +2mp (2t)) 

(37) 

= {K z (j-j' + 2m P ;t) + K z (j+j'-l + 2mp;t)), 



m= — oo 



in agreement with (34). 

The trace can be constructed, but, again, yields nothing new and so I will not 
write it out. 

9. The D-N interval 

Having treated the D-D and the N-N intervals, I now complete the set by tack- 
ling the hybrid, D-N case, although nothing especially new is expected to emerge. 

I choose j = as the D end point, i.e. y(0) = 0, and set the mirrors at 
j = and at j = p + 1/2 enclosing p free vertices, as in the N-N case. The 
images are located at j — 2m(2p + 1) for even reflections and at the opposite values, 
—j — 2m{2p + 1), for odd, (— oo < m < oo). The even terms enter with the sign, 
(— l) m , while the odd ones have (— l) m+1 . Hence the heat-kernel is, 

K» D (j,j';t)= (-l) m U z (j,j' - m (2p+l);t)-K z (j,-j' -m(2p+l);t)\ 

m= — oo ^ ' 

= E (- 1 ) m (^z0' - f+ m ( 2 P + !); *) ~ K zU + f+ m ( 2 P + !); *)) • 

m= — oo ^ ' 

= e E ( _ Ij-j'+m(2p+l)(2t) — Ij + j> +m (2p+l)(2t) J . 

m= — oo ^ ' 

(38) 

The transformation of this into the eigenfunction form provides a use for the 
particular twisted identity, (16), corresponding to anti-periodic functions. I repeat 
it here, with p — )■ 2p + 1, 

oo 2p 
p~ z ( -\\ m T (A- -2zsin 2 hr(2n+l)/2(2p+l) 

e 2^ l > i i+m(2 P +i)^J - „ 1 2^ e ' 

m=-oo n=0 (39) 

7r(2n + 1) j 

X COS — ; — . 

(2p+l) 
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Substitution into (38) produces, at first, 



1 



2p+l 



p-1 2p 
n=0 n=p+l 



-4tsin 2 (7r(2n+l)/2(2p+l)) 



os ,(2n + l)(,-/)_ cos ,(2n + l)( J+/) 



2p + l 



2p+l 



(40) 

after noting that the n = p summand is zero. As earlier, the two sums are equal 
and so, 



A P_1 / \ 

K ND (j i'-t) = V^ e -4tsin 2 (7r(2n+l)/2(2p+l)j 

V \JiJ i ) 2 p _|_ 1 



n=0 



. n(2n + l)j . n(2n + l)j' 

x sin — — sin — — , 

2p + l 



(41) 



2p + 1 

from which the eigenvalues and eigenfunctions can be read off, in a standard way, 
as, 

Afr J -4sm 2 ^ tv-, n = 0, ...p-l 



2(2p+l) 



and 



v5 N U) 



2 7r(2n + l)j 

sin 



V2p~fT ' 2p + l 
These agree with available expressions e.g. [4], [6]. 



10. Lattice Green functions 



In the continuum, the Laplace transform of the heat-kernel is the Green func- 
tion, and this holds in the discrete case too. 

The transform of the right-hand side of (15) (times 1/2) is the Green function 
for the p-cycle, C, 



p-i 



1 p 2ni(n+a)(j-j')/p 

G, C (i./.«.7) = j-E- 



4p ~ sin 2 (77(71 + a) jp) + sinh 2 7 
= G p (j -/,0,a,7) 
= G p (j ~ 3',a,l) ■ 



(42) 



In order to apply the transform to the left-hand side of (15) the sum has to be 
arranged so that the order of the Bessel function is always non-negative i.e. m is 
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such that j + mp > 0. Those terms for which j + mp is negative can be converted 
to positive j + mp using I_ j _ mp = I j+mp . 

It is best to use residues mod p and set j = Jp + r with J e Z and < r < p. 
The positive terms then correspond to m > — J + 1 and the negative ones to 
m < — J — 1. The zero term, m = — J, gives an I r and the sum therefore can be 
rearranged to, 

oo 

1 Jp+r+mp T~ e 1 Jp—r+mp) ~T~ e J rl z J ) 

m=-J+l 

allowing the Laplace transform to be taken giving, 
1 

2 sinh 27 

'm=-J+l 

(43) 

The m-summation can be done as before and, taking the real part, yields the 
identity, 

p— 1 

1 cos2n(n + a)j /p cos 2ir oJ f cosh 2^r sinh 2p^ 
— / — 9~7 \ 9 — = ; ; sirm 27 r 

sin27ro;J sin 27ra sinh 2^r 



^ \ —2-fmp ( p—2-fJp ( —2-fr —2irimct , 2^r 2-rrima+4:TTiJa\ , 2nija —2-fr 



sinh 27 cosh 2p^ — cos 27rct 
cos(27rctJ) t/ p _ 1 _ r (cosh27) + cos (2tv a(J + 1)) U r _ 1 (cosh2-f) 
2 (T (cosh 27) - cos 2na) 

(44) 

again in terms of the Chebyshev polynomials, which are the basic building blocks 
for the Green functions of the free equation. 

As a check, when j = (J = = r) equating this to the real part of (42) 
gives the identity, (18). Also, when a = 0, it is periodic in J of period 1 (i.e. in 
j of period p) as it should be. Actually it is then independent of J the expression 
being 12 

^ , \ cosh ( 2r — 0)7 
Re G c JJp + r, 0,7)- 



2 sinh 27 sinhj>7 
U p _ r _ 1 (cosh2 1 ) + U r _ x (cosh 2 7 ) ( 45 ) 

2(T p (cosh2 7 )-l) ' 



12 



This is the standard result found in Feller, [11]. See Montroll and Weiss, [33], equn.II.19. 



15 



I also write out the anti-periodic Green function, (a = 1/2), 



= J ^_ r _ 1 (cosh2 7 )-^_ 1 (cosh2 7 ) ( 46 ) 
2(T p (cosh2 7 ) + l) 

When a / 0, the twisted periodicity appears in its SO (2) aspect and, as a 
by-product, yields the imaginary part of (42) which can be confimed by performing 
the summation, (43). I will not expose the easily found expressions. 

This lattice Green function has also been evaluated using Laplace transforms 
by Cojocaru, [20], equn.(20). 

The formula (44) generalises that of Bendito, Encinas and Carmona, [34], who 
give the untwisted (45). 

The interval Green functions, which are real, can be obtained by combining 
the cycle expressions, (45) and (46), according to the image forms, (29), (37) and 
(38). 

The Dirichlet Green function for the p-path, V, is found to be, after some 
simple trigonometric rearrangement, 13 

D( , , _ ^- r (cosh2 7 )^_ 1 (cosh2 7 ) 

pl ' ' 7J t/- p (cosh2 7 ) 1 ] 

where r > r' . If r < r', the symmetry of the Green function in r and r' can 
be employed. This result agrees with that of Chung and Yau, [1], derived in a 
more involved way, and that of Bendito, Encinas and Carmona, [34], obtained 
via a Sturm-Liouville approach which illuminates the structure of (47) and avoids 
mention of the eigenproblem. 

The Neumann Green function is found to be, again after simple trigonometric 
manipulation, 

nN( , s ^- r (cosh2 7 ) y r ,_ 1 (cosh2 7 ) 

^ (r - r - 7)= 2(cosh2 7 -l)^_ l( cosh2 7 ) - (48) 

which agrees with [34]. The third kind Chebyshev polynomial, V, arises through 
the Neumann offset r — > r — 1/2, r' — > r' — 1/2. 



13 On the interval j = r, j' = r' . As a tactical point, I find it easier to use trigonometric formulae 
than relations between Chebyshev polynomials. 
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Finally I find the N-D Green function to be, 



G nd , r , , _ V r (cosh2 7 )^_ 1 (cosh2 7 ) 

^ ^ r ' 7j " F p (cosh 2 7 ) ' " r • [ yj 

To attain this form, one sees from (38) that it is the anti-periodic cycle Green 
function, (46), that is required here. 

Bass, [35], gives the Green functions for all three interval types, but has a more 
complicated structure. He obtains the form for the Dirichlet and Neumann cases in 
terms of the periodic one by algebra. 



11. Graph considerations 



The stepping operator E is the coordinate (or vertex) space representation of 
the abstract stepping operator E, and E + E _1 is the adjacency operator, A, having 
the usual adjacency matrix representation on the integer lattice, Z, coordinates. 

One sees from the definition (4) that K(£) is a generating function for the 
powers of the adjacency operator, 

K Z (*) = , 

with 

K z (jj';t) = (j\K z (t)\j'). 
A similar relation holds on the p-cycle projection, C p = Z/pZ. According to 

(6), 

K c (t) = e A c\ (50) 

which are operators in the p-dimensional vector space describing C p . A c is repre- 
sented by the usual periodic adjacency matrix. 
The trace of the expansion of (50) yields, 

00 , 

trK c (t) = ]T-trA£, (51) 

k=0 

and, from the definition of adjacency matrix, tr is the number of closed paths 
(circuits) of length k on the p-cycle, which equals, by homogeneity (or transitivity) 
p times the number of circuits starting and finishing on any one vertex. Denoting 
this latter number by g k , one can derive identities from the various forms of the 
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left-hand side of (51) which allow the computation of g k , or, equivalently, of the 
series generating function defined by, 

oo 
k=0 



The Laplace transform of (51) produces, 

1 f°° 

g(a) = - / dte-^tTK c (t), 
& Jo 



which is the trace of the cycle lattice Green function, already computed, in (18) or 
(45), the relation between variables being 2 cosh 27 = 1/a > 2. Hence one has the 
explicit rational formula for the generating function, 

= 1 sinh2p 7 1 = 1 flp_i(l/2(T) 

9[a) 2a sinh27 cosh2p7 - 1 2a (T p (l/2a) - l) ' 

A more general quantity in graph theory is the number of circuits of length k 
with / bumps (or immediate backtrackings), attached to a vertex. Denoting this 
number by f kh the corresponding generating function is, 

00 

f(u,a) = fikU l a k , 

k,l=0 

so that, as a special case, g(a) = /(l,cr) and for the number of circuits with no 
backtrackings, f(a) = /(0,<r). 

If the graph is <i-regular (for the cycle d = 2) there is the relation, [36], 

f{1 - U - a) - 3 ( , , ' ■ n ), (53) 



1 — u 2 a 2 \ 1 + u(d — u)a 2 ) 

which allows the more interesting f(u,a) to be obtained from the simpler g(a). 
Trigonometric expansion (e.g. Loney, [22] §§293-294, can be used to evaluate this 
but it is simpler to employ Chebyshev polynomials as rapid computation is available 
through, 

U n (x) = tv (C n (x) Q) , T n {x)= l -trC n {x), V n {x) = tr (C n (x) R) , 



where C, Q, R and S are given by, 



14 



,0 1\ ^ /0 0\ /0 0\ n ( -\ 1 

CV) =I -1 2, • 1 • R= -1 1 ' S= -1 1 



14 I introduce S 1 because the denominator of the Neumann Green function, (48), can be written 
U p + U p _ 2 -2U p _ 1 =tr(CPS). 
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Hand evaluation quickly produces the first two examples, 

l + (l-u)a 1 + (1 - u 2 )a 2 

hM= l-{l + u)a > f * M= l-(l + U )W 

where the subscript refers to the number of vertices, i.e. p. 

For larger p, it is sufficient to give the g p from which the f p (u, cr) can readily 
be found on applying (53). I list the results for p = 1 to p = 12, 

1 1 cr - 1 2 a 2 -1 a 2 + a- l 



2(7-1' 4(7 2 -l'2a 2 + a-l , 4a 2 -l' 2 a 3 - 3 a 2 - a + 1 ' 

3(7 2 ~1 (7 3 ~2(7 2 -(7 + l 2(7 4 - 4(7 2 + l 

4 (7 4 - 5 (7 2 + 1 ' 2(7 4 + 3(7 3 - 4cr 2 -(7 + l' 8(7 4 - 6(7 2 + l' 

a 4 + 2a 3 -3(7 2 -(7 + l 5a 4 -5(7 2 + l 



(54) 



2(7 5 - 5(7 4 - 4(7 3 + 5(7 2 + (7 -l , 4(7 6 - 13 (7 4 + 7(7 2 - 1 ' 

The ensuing f p (u,a) agree with Bartholdi, [36], (who uses more involved alge- 
bra) except that his even expressions appear to have misprints. 

Taylor expansion allows the coefficients f lk to be read off and the correctness 
of f 4 (u,a), say, can be checked. For interest, I write it out, 

f 4 (u, a) = 1 + 2 u a 2 + 2 (u 3 + u 2 + u + l) a 4 

+ 2(m 5 +2« 4 +4« 3 + 6« 2 + 6«) (t 6 + .... 

It is possible to confirm numerically the determinantal identity (17) by comput- 
ing the Laplacian determinant using the explicit forms of the generating functions 
(54). This is no more than a check of algebra of course, but it makes a comforting 
exercise. 

The Laplacian graph eigenvalues, denoted A n , (n = 1 — >■ p), are related to the 
adjacency eigenvalues, \i n by A n = 2 — fj, n , in this cycle case. The first A n , i.e. X 1 , 
is zero and so fi x = 2. The log Laplacian determinant, defined omitting the zero 
eigenvalue, is, therefore, 

v v 
logdet£(p) = ^logA n = ^log(2-^J 



n=2 n=2 

oo 



(p-l)log2-££--^ 



k=l n=2 



(p-l)log2-£^(trA^-2 fc ) 



(55) 



k=i 

oo 



(p-l)log2-£-l r ( P ( 7fc -2 fc ). 
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Now define the subtracted generating function, g*, by 

which is finite at a = 1/2, then (55) becomes, 

det' L (p) = 2^ 1 exp ( - ^ (p 9 ;(a) -p + l)) (56) 

which can be evaluated and gives the right answer 15 , det^(p) = p 2 , case by case, 
using (54). The integration is the reverse of the step taking (17) into (18) for 
arbitrary p. 

A further check of the expressions, if one were needed, is provided by computing 
the graph adjacency matrix from the expansion of the Green function in powers of 
p = cr/(l — 2a). The coefficient of — p 2 is the combinatorial Laplacian matrix, L, 
from which the adjacency matrix can be found. The definition of Laplacian I use is 
the operator relation, G = p/il + pL). 

I do not give the details since the answers are as expected. However, for the 
p-path, I do write out the Laplacian matrices for the D-D, N-N and D-N boundary 
conditions computed from (47) , (48) and (49) , using p = 3 as an illustration, 

/ 2 -1 \ / 1 -1 \ / 2 -1 \ 

L D =l-l 2 -1 lL N = -1 2 -1 lL DN = -1 2 -1 . 
V -1 2 / \ -1 1 / \ —1 1 / 

(57) 

The adjacency matrix is obtained by replacing the diagonal elements by zeros. 

The Neumann diagonal elements correspond directly to the degrees of the free 
vertices whereas the first and last Dirichlet elements include the effects of the two 
fixed end points, which are not part of the free graph and serve only to make L D 
non-singular. 16 This is confirmed by the D-N matrix. 

A combinatorial expression for the number of closed circuits, g k , on the p-cycle 

15 According to one version of Kirchhoff's matrix-tree theorem, it gives the number of spanning 
trees, det£(p)/p, correctly as p for the cycle. 

16 In the terminology of Chung and Yau, [37], the graph is an induced subgraph of a host graph, 
which includes the end, or boundary, points. This can be seen from the form of the matrices in 
(57). L D is obtained from L N , in two higher dimensions, by striking out the rows and columns 
that refer to the end points. 
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also follows directly from the operator definition, 




m — — [k I p] 
k — mp even 



P 



1 



(k 



trtE + E-y, E p = I 



mp)/2 



k 



) 



12. The thermodynamical limit 

As mentioned before, the covering space (an infinite lattice) corresponds to the 
thermodynamic limit of the covered (a finite lattice) (e.g. [4]). For example, very 
easily, the Laplace transform of the heat-kernel on Z, (5), gives the thermodynamic 
limit of the cycle free Green function on Z/pZ directly as, 



The parameter, 7, can be extended into the complex plane so long as it retains a 
positive real part (assuming j > j') and then one can compare with the formulae in 
Katsura and Inawashiro, [38], Appendix, (although I have problems with equation 



13. Remarks 

Most of the preceding expressions and relations are rather elementary examples 
of more general structures that have been investigated for many years, in many 
different areas, as mentioned in the introduction. For example, (bulk) lattice Green 
functions, in all dimensions, have been rewritten in terms of Bessel functions both 
for numerical and formal reasons, an early reference being van der Pol and Bremmer, 
[39], p. 368, in connection with electrical networks. Later works include [40] and [38] 
where Mellin-Barnes type integrals for Bessel products are utilised to give more 
explicit representations of the lattice Green functions, a topic of some considerable 
activity and manipulative complexity. See the reviews, Guttmann, [41], and Zucker, 



Incidentally, van der Pol and Bremmer also take the infinite square lattice 
expression to an elliptic integral, see [39], equn.(91). A number of forms for this 



G(j, j\ 7) 



e -2 7 (j-j') 




2sinh27 ' 



(A.4)). 



[42]. 
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quantity have been obtained over the years, the earliest being, perhaps, by McCrea 
and Whipple, [43] , for a random walk problem. More mathematical are the compu- 
tations by Stohr, [44], who touches on the resistor network realisation which itself 
has also seen a goodly amount of development. Flanders, [45], lists some square 
expressions. 

14. Potentials, higher dimensions and continuum limit 

I wish to mention possible extensions to the preceding calculations which I will 
not be carrying out here. 

The analysis in this work has all been concerned with free propagation. An 
obvious development would be to consider the case with a potential described by 
the recurrence, 



instead of (1). In order to project down to Z/pZ, for example, V(j) would have to 
be periodic, V(j + p) = V(j). Working on the interval, V, would also be possible 
directly and related to a Sturm-Liouville problem. 

Unless exact solutions were available, one might be interested in the discrete 
analogue of the short-time behaviour of the heat-kernel, and the resulting expansion 
coefficients. One possible structure has been analysed by Iliev, [46]. 

I have also restricted myself to one dimension. Higher (free) dimensions can 
be assembled in the usual way by taking the cartesian products of one-dimensional 
structures, for example a torus from the circle, a rectangle from the interval and and 
cubic lattices from the integers. Some basic notions are outlined in [1] and a more 
extensive discussion of tori in Chinta, Jorgenson and Karlsson, [47], [48] (see also 
Chaumard, [49]), who are particularly concerned with the continuum limit, another 
important topic I have not addressed. 

15. Conclusion 

I have made a slight extension to the computation of the free heat-kernel on 
a p-cycle and have used the resulting flexibility to derive a bigger variety of Bessel 
relations, for example. 

An organisational point is that, if one is interested just in the discrete resolvent 
(the left-hand side of (18)) then an appeal can be made to the identity (17) and 
the Bessel form can be bypassed. 




(58) 
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I have also calculated the heat-kernel on a discrete interval for various boundary 
conditions using both modes and images. As usual, the Neumann images are all 
positive while the Dirichlet ones alternate in sign (positive for even reflections and 
negative for odd). In the continuum, when adding Neumann to Dirichlet the odd D 
images cancel corresponding N ones leaving just the images for a periodic circle of 
circumference twice the size of the interval, as is well known. This cancellation does 
not happen in the discrete case due to the offset of 1/2 to the Neumann mirrors. 

Expressions for the free Green functions in terms of Chebyshev polynomials 
are derived. That on the D-N path is perhaps novel. As a check, the combinatorial 
Laplacian matrices are extracted. 

Relatedly, rational generating functions for circuits of specific length and num- 
ber of bumps on the cycle graph are readily computed and compared with those of 
Bartholdi, [36]. 

Paths and cycles are graphs and these calculations illustrate very general 
graph-theoretic procedures. I intend to expand on this aspect elsewhere. 
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